Beliaev damping of quasi-particles in a Bose-Einstein condensate 
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We report a measurement of the suppression of collisions of quasi-particles with ground state 
atoms within a Bose-Einstein condensate at low momentum. These collisions correspond to Beliaev 
damping of the excitations, in the previously unexplored regime of the continuous quasi-particle 
energy spectrum. We use a hydrodynamic simulation of the expansion dynamics, with the Beliaev 
damping cross-section, in order to confirm the assumptions of our analysis. 



In a Bose-Einstein condensate (BEC) collisions with 
distinguishable excitations (impurities) and collisions 
with indistinguishable excitations (quasi-particles) differ 
profoundly. These differences follow from the quantum 
exchange symmetry of the indistinguishable excitations, 
and from a different excitation spectrum. 

Impurity collisions within a BEC have been measured 
previously. Using Raman spectroscopy the microscopic 
onset of superfluidity was measured and found to be in 
general agreement with prediction W. The collisional 
dynamics and interaction between two distinguishable 
slowly moving condensates was measured pi and found to 
agree with simulation H . Macroscopic superfluid behav- 
ior has also been demonstrated, involving the interaction 
of the condensate with large-scale optical dipole-potential 
structures ||], (§• 

The case of identical particle collisions has been exten- 
sively studied using many-body theory, starting with || . 
Recently, these results have been applied to BEC explic- 
itly @ > H] i |§ ; (lO) • I 11 this letter we present a measure- 
ment of collisions between quasi-particles and the BEC, 
at velocities near and above the superfluid critical veloc- 
ity v c . 

According to the Fermi golden rule, the rate of scat- 
tering within a homogenous BEC is given by [0|: 
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where a is the s-wave scattering length. The wavenum- 
ber k of the excitation is in units of = V8 irna, the 
inverse healing length of the condensate. The free parti- 
cle velocity of the excitations is v k = hk^ 1 /m, where m 
the mass of the BEC atoms, n is the density of the con- 
densate. The superfluid critical velocity is v c — y/ n/m, 
where fj, = gn is the chemical potential of the BEC, and g 
is A-Kh 2 a/m. The integral is over all possible momentum 
transfers q, and all angles Q. The <5-function requires en- 
ergy conservation between the initial energy Ei and the 
final energy Ef after collision, where we express energy in 
units of /i. The factor A q -k is the q-dependent, momen- 
tum conserving, scattering matrix element of an excita- 
tion of wavenumber k, which may include suppression or 
enhancement of the collision process due to many-body 
effects. 

For impurity scattering within a condensate [ pT[ , 
|A g; fe| 2 is given by the structure factor S q [Q. The ini- 



tial energy Ei — E%, where E® = k 2 is the impurity 
dispersion relation, for impurities with mass m. The fi- 
nal energy for this process is Ef = E± + ^ q , where 
w v , , 2q 2 is the recently measured Bogoli 



ubov dispersion relation |14| . 

Below v c collisions between the impurity and conden- 
sate are completely suppressed (see Fig. I) by the 5- 
function requiring conservation of energy and momen- 
tum. 
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FIG. 1. Cross-section for collisions in a homogenous con- 
densate. The momentum is in units of the inverse healing 
length, £ = V8irna. The cross-section is in units of the free 
particle scattering cross-section for identical particles 87ra 2 . 
The solid line is the theoretical curve, from Eq. (3) for 
quasi-particles travelling through the condensate. The dashed 
line is the cross-section for impurity scattering within the con- 
densate [1] , which vanishes below the superfluid critical veloc- 
ity v c . The inset shows the allowed momentum manifold for 
quasi-particle collisions due to conservation of energy and mo- 
mentum. The gy and q± axes correspond to the parallel and 
orthogonal components of the scattered momentum respec- 
tively, where tan(8) = q±/q\\. The manifolds represent the 
experimental fc's 4.09 (outermost), 2.63, 1.73 and 1.06 (inner- 
most). 

In the case of quasi-particle collisions only the Bogoli- 
ubov dispersion relation is relevant, and the S- function 
in Eq. (@) is simply 5(E% - E^ - £f_ q ). We solve this 
condition, and find the angle 6 between the initial direc- 
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tion k and the scattered direction q, to be (see inset of 
Fig. 1): 
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This result differs in a qualitative way from impurity 
scattering, since Eq. (|J) has solutions for any finite k. 
There is no longer any well-defined critical velocity at 
which collisions are completely suppressed. However, not 
all angles are allowed. At a given k we find that the 
maximal allowed angle is cos(9 max ) = \J (fc 2 + 2)/2/(fc 2 + 
1). At the limit of small fc, this angle approaches zero, 
and collisions are allowed only for q parallel to k. 

The appropriate suppression term |A 9; fc| 2 for quasi- 
particles has been calculated Q], f|. In this work we 
expect mainly Beliaev processes which involve creation 
of lower energy excitations. The Landau damping rate 
is expected to be an order of magnitude slower than the 
observed Beliaev collision process || . 

We start with the atomic interaction Hamiltonian 
H' = ^ Y,j,i, m ,n a j a l a ma>n8j+l-m-n, where V is the 
volume of the BEC, a\ and a, are the atomic creation 
and annihilation operators at wavenumber i. We approx- 
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the condensate. We take the Bogoliubov transform al = 
(u p bp — Vpb^p), with u p and v p the appropriate quasi- 
particle amplitudes, which were recently measured Jl5| . 
We are interested in terms of the form bkb\._^> q , that re- 
move a quasi-particle of wave number fc, and create two in 
its stead. Calculating the matrix element prefactor of this 
term in the atomic interaction Hamiltonian, we arrive at 
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This result can be viewed as the explicit zero tempera- 
ture limit of more general calculations [fLOl. 

Applying Eq. (H) and |v4 g; fc| 2 to Eq. <^), and using the 
Feynman relation Jl^] S q - .._, ( 
of excitation-condensate collisions: 



E®/Ey, we arrive at the rate 
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The effective cross-section for the quasi-particles, 
is shown in Fig. 1 (solid line). For large fc, <jjf ap- 
proaches 87ra 2 , compared to Ana 2 for impurities (dashed 
line). This enhancement by a factor of 2 is due to the bo- 
son quantum mechanical exchange term. In Eq. (H), for 
small fc, we verify that the scattering rate indeed scales 
as fc 5 , which is the classic result Hj, ||. In particular, it 
remains finite even for v k < v c , in contrast with impurity 
scattering. 

In M, the identical particle collision cross section for 
large fc was measured to be 2.1(±0.3) x Ana 2 . Scattering 



rates in four-wave mixing experiments in BEC |17|] were 
also shown to agree with the high-fc limit of Eq. ([|) [Q . 

In the opposite regime of extremely low wavenumber, 
where the energy levels are discrete, Beliaev damping 
was observed for the scissors mode of a BEC Q. The 
discrete energy levels were tuned so that Beliaev damping 
of the initial mode to exactly one mode of half the energy 
was achieved. Eq. (||) did not apply, since there was no 
need to integrate over various scattering modes. 

Our experimental apparatus is described in [ jl3| . 
Briefly, a nearly pure (> 95%) BEC of 10 5 s7 Rb atoms 
in the \F, m/) = |2, 2) ground state, is formed in a QUIC 
type magnetic trap |20[| . The trap is cylindrically sym- 
metric, with radial (f) and axial (z) trapping frequen- 
cies of 2ir x 220 Hz and 2tt x 25 Hz, respectively. Thus 
£ = 0.24/im via averaging in the local density approxi- 
mation (LDA) p[. 
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FIG. 2. Absorption TOF images of excited Bose-Einstein 
condensates, (a) Absorption image for k = 2.63, with the 
large cloud at the origin corresponding to the unperturbed 
BEC. A clear halo of scattered atoms is visible between the 
BEC and the cloud of unscattered outcoupled excitations, (b) 
Absorption image for k = 1.06. For this value of k the dis- 
tinction between scattered and unscattered excitations is not 
clear, since both types of excitations occupy the same region 
in space. 

We excite quasi-particles at a well-defined wavenum- 
ber using two-photon Bragg transitions The two 
Bragg beams are detuned 6.5 GHz from the 5<Si/2, F = 
2 — ► 5P3/2,F' = 3 transition. The frequency difference 
Aw between the two lasers is controlled via two acousto- 
optical modulators. Bragg pulses of 1 msec duration are 
applied to the condensate. The angle between the beams 
is varied to produce excitations of various fc along the 
z-axis, thus preserving the cylindrical symmetry of the 
unperturbed BEC. The beam intensities are chosen to 
excite no more than 20% of the total number of atoms in 
the condensate. 

After the Bragg pulse, the magnetic trap is rapidly 
turned off, and after a short acceleration period the in- 
teraction energy between the atoms is converted into 
ballistic kinetic energy j22|. After 38 msec of time-of- 
flight (TOF) expansion the atomic cloud is imaged by 
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an on-resonance absorption beam, perpendicular to the 
z-axis. Fig. 2a shows the resulting absorption image for 
k = 2.63, with the large cloud at the origin corresponding 
to the BEC. A halo of scattered atoms is visible between 
the BEC and the cloud of unscattered outcoupled excita- 
tions. No excitations with energy greater than that of the 
unscattered excitations are observed, confirming our low 
estimate of the Landau damping rate. Fig. 2b shows the 
absorption image for k — 1.06. For this k value the dis- 
tinction between scattered and unscattered excitations is 
not clear in the image, since both types of excitations 
occupy the same region in space. 

At a given k the number of excitations is varied by 
scanning Aui around the resonance frequency oJu ■ The 
number of excitations N mom is measured by determining 
the total momentum (in units of the recoil momentum 
hk^ 1 ) contained in the outcoupled region outside the 
unperturbed BEC. This region includes all the scattered 
and unscattered excitations, in the direction of k. Ther- 
mal effects are removed by subtracting the result of an 
identical analysis over the other side. The results are 
shown in Fig. 3a. 
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FIG. 3. Quantifying the amount of collisions measured for 
k = 2.63. (a) Measured momentum of the outcoupled atoms 
vs. Aui. The momentum is in units of the recoil momentum, 
and is normalized by the total number of atoms, (b) The 
measured ratio between N mom and Ncount (both defined in 
the text). Every collision outcouples more atoms, increasing 
N coun t but leaving N mom almost unchanged. 

In order to quantify the amount of collisions, despite 
the lack of separation between scattered and unscattered 
excitations, we take the ratio between N mom and the 
counted number of atoms N count , in the same region. 
The resulting N mom /N count , as a function of Aw, are 
shown in Fig. 3b. The ratio is seen to be independent of 
the number of excitations, and appears to be an intrinsic 
property of a single excitation. Each collision between 
an excitation and the condensate creates an additional 
excitation that is counted in the outcoupled region, in- 
creasing N count , while the momentum (N mom ) in this 
interaction is conserved. Thus the ratio, N mom /N count , 
is a good quantifier of the amount of the collisions, even 
at low k 1 23 1 . Bosonic amplification of the collision rate 
jyj , which would appear as minima in Fig 3b, is not 
observed. 



At a given k we define the overall probability for an ex- 
citation to undergo the first collision If we ignore sec- 
ondary collisions the result is N mom /N count = 1/ (l+Pk), 
since each collision outcouples, after TOF, an additional 
particle pifl . Using this relation we infer p e ^ p for the 
various measured fc's. 

We expect the scattering probability pk to be equal to 
Vkteff, where t e ff is the effective interaction time 
of the excitation with condensate and h is the average 
density. 

We assume t e ff to be fc-independent, divide pk by Vk, 
and arrive at a value that is proportional to the scat- 
tering cross section (since the h is constant for all k). 
This assumption will be tested below, but must be valid 
for sufficiently low k 7 for which the TOF expansion low- 
ers the density rapidly, turning off collisions before the 
excitations move significantly. 

The ratio, p e ^ p / '(hvkt e f /) is shown in Fig. 4 (■) and 
seen to agree with the theoretical suppression (solid line) 
calculated as a LDA average [|l2| of Eq. (||) |25| . Since 
t e ff is not known, absolute calibration is not possible. 
Therefore, the results are shown with arbitrary units. 

In order to verify the validity of these simplifying as- 
sumptions, 3-D simulations are performed, using the hy- 
drodynamic Gross-Pitaevskii equations pfij ]. The simu- 
lations [^2| include the TOF dynamics. The expansion 
of the main BEC cloud is taken from expressions in [^7) , 
computed for an elongated condensate by the hydrody- 
namic equations. The excitations travel within the time 
varying mean field potential created by the expanding un- 
perturbed condensate. The excitations collide with the 
expanding condensate with the correct local scattering 
cross section, including the angular dependence, taken 
from Eq. ([|). The resulting distribution of simulated 
outcoupled atoms is analyzed by the same method as the 
experimental data. 

We plot the simulated Pk/vk (dashed line, in the same 
units as experiment). The simulation and Beliaev damp- 
ing theory agree in the regime of experimentation, in- 
dicating a window of validity for the assumptions used 
in analyzing the experimental data. At large k (above 
k = 4.09), t e ff is made shorter by the rapid transit of 
the condensate by the excitations. At low k (less than 
k = 1.06) many of the collisional products remain inside 
the condensate volume, and are not counted, preventing 
analysis in this regime. 

The arbitrary unit suppression used in Fig. 4 (both 
for simulation and experiment) was set at k — 2.63. The 
absolute cross-sections obtained by comparing the ex- 
perimental data to the hydrodynamical simulations were 
higher than the theoretical values by a k-independent 
overall factor of 2.36 ± 0.08, which is not understood. 
This factor may be caused by various inaccuracies in the 
TOF parameters of the simulation. However, the trend in 
the experimental analysis is robust, and does not depend 
on absolute calibration. 
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FIG. 4. Suppression of identical particle collisions. Scat- 
tering cross-section is in arbitrary units. The momentum is 
in units of the inverse healing length after LDA averaging, 
£ = 0.24/im. The theoretical curve (solid line) is a LDA aver- 
age of Eq. (3) [25] . The assumptions of our analysis are tested 
using hydrodynamic simulations, and found to agree with 
Beliaev damping theory in the experimental regime (dashed 
line) . 



We also set the collision rate artificially to zero in the 
simulation, and find N mom /N count to be unity within 2%, 
for all k. This implies that there are no significant other 
mean- field repulsion effects along the z-axis p2]| , confirm- 
ing our assumption of momentum conservation. 

In conclusion, we report a measurement of the suppres- 
sion of the collision cross-section for identical particles 
within a Bose-Einstein condensate. We find the suppres- 
sions in our experiment in agreement with a calculation 
of Beliaev damping rates, within an overall factor. We 
use a hydrodynamic simulation of the expansion dynam- 
ics, in order to verify our analysis of the experiment. 
This represents the first measurement of this effect in 
the quasi-particle continuous spectrum regime. 

In the low temperature limit discussed in this Let- 
ter, the non- linear coupling mechanism of the interac- 
tion Hamiltonian is predicted to generate squeezing and 
entanglement of quasi-particle excitations pq | . Further 
experimental study should also lead to an understanding 
of these processes in the context of the decoherence of 
excitations and coherent states due to collisional losses 
within a BEC ||, @, fl). 
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